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ABSTRACT 



Context. Stellar parameter determination is a challenging task when dealing with galactic giant stars. The combination of dif- 
ferent investigation techniques has proven to be a promising approach. We analyse archive spectra obtained with the Short- 
Wavelength-Spectrometer (SWS) onboard of ISO, and new interferometric observations from the Very Large Telescope MID -infrared 
Interferometric instrument (VLTI/MIDI) of a very well studied carbon-rich giant: TX Psc. 

Aims. The aim of this work is to determine stellar parameters using spectroscopy and interferometry. The observations are used to 
constrain the model atmosphere, and eventually the stellar evolutionary model in the region where the tracks map the beginning of 
the carbon star sequence. 

Methods. Two different approaches are used to determine stellar parameters: (i) the 'classic' interferometric approach where the effec- 
tive temperature is fixed by using the angular diameter in the A'-band (from interferometry) and the apparent bolometric magnitude; 
(ii) parameters are obtained by fitting a grid of state-of-the-art hydrostatic models to spectroscopic and interferometric observations. 
Results. We find a good agreement between the parameters of the two methods. The effective temperature and luminosity clearly 
place TX Psc in the carbon-rich AGB star domain in the H-R-diagram. Current evolutionary tracks suggest that TX Psc became a 
C-star just recently, which means that the star is still in a 'quiet' phase compared to the subsequent strong-wind regime. This is in 
agreement with the C/O ratio being only slightly larger than 1 . 



Key words. Stars: AGB and post-AGB 
metric - Techniques: spectroscopic 



1. Introduction 



Stars: atmospheres - Stars: carbon - Stars: fundamental parameters - Techniques: interfero- 



The asymptotic giant branch (AGB) is the late evolutionary stage 
of low- to intermediate- mass stars (1 - 8 Mq). On the early - 
AGB the carbon-to-oxygen-ratio is smaller than one. After sev- 
eral thermal pulses, the atmospheres of objects with masses be- 
tween 1 - AMq will very likely turn from o xygen-rich into 
carbon-rich because of the third dredge-up (llben & Renzinil 
Il983l) . Good estimates of stellar parameters are needed for a 
profound understanding of the evolution of this stage. Their 
determination is a challenging task because of the complex- 
ity of the atmospheres of these objects. It is demonstrated that 
the combined use of spectroscopic and interferometric observ- 
ing techni ques can efficiently help to ascertain stellar parame- 
ters ( e.g. IWittkowski et aDl200ll |2008[ 1201 \\ [Neilson & LesteH 
2008; Pala dini et al.l201 ltlSacuto et alJ201 1 at Marti- Vidal et all 



Catalogue of Variable stars (Sa mus et al .1 120091) with a mean 
brightness of ~ 5 mag and a peak-t o-peak amplitude of 0.4 mag 
in the V-band dJorissen et al.ll20lTI). Dis t ance estimates range 
from 275 to 315 pc (Ivan Leeuwenll20"07t IClaussen et al.l fl987: 



Bergeat & Chevallier 



2011). At the same time these observations provi de constraints 
for ex i sting model atmosp heres: e.g. C OM ARCS (| Aringer et al.l 



120091) . iHofner et al l d2003h . PHOENIX (Hausc hildt et alJ 1999), 
ATLAS (e.g. lKurucz)1993l) . CODEX (llreland et alJ2008Ll2011h . 



TXPsc is one of the brightest and closest carbon-rich AGB 
stars. It is listed as an irregular variable in the General 



* Based on observations made with ESO telescopes at Paranal 
Observatory under program IDs 74.D-0601, 60.A-9224, 77.C-0440, 
60.A-9006, 78.D-0112, 84.D-0805 



2005TT Different ISO/SWS spec tra of 
TXPsc dJ0rgensen et al.l2000l;lGautschv-Loidl et al.l2004 show 
that there is a difference in the 3^m feature between 1996 
and 1997. With the he lp of plane-parallel, hydrostatic models 
J0rg ensen et al.l (l2000h interpret this difference as a change in 
temperature of ~100 K. The photometry and spectr a were suc- 
cessfully modeled by Gautschv- Loidl et al.l d2004l) with dust- 
free dynamical models that reproduce the region between 1 - 
5 /im. They claim that the region between 8-9 //m is affected 
almost solely by CS. The authors suggest that observations from 
8-9 fim show the deep photosphere, while the expected absorp- 
tion of HCN and C2H2 originating from the higher layers is not 
observed. As this object is almost (carbon) dust-free we may in- 
fer it became a carbon star quite recently. Therefore, this star is 
a perfect candidate to constrain the region where the transition 
from oxygen- to carbon-rich occurs. 

In this work we present a study of the atmosphere and a stel- 
lar parameter determination for TX Psc. We combine spectro- 
interferometric observations of VLTI/MIDI and spectroscopic 
observations from ISO/SWS and compare them to geometric, 
hydrostatic, and evolutionary models. 



1 



Klotz et al.: Constraining stellar parameters for TXPsc 



100 



-50 





i i i i i i i i i 


- \ 


- 

























-100 -50 50 100 

u [orcsec"'l 

Fig. 1. iV-band spectrally dispersed Mv-coverage of the MIDI ob- 
servations of TXPsc. Colour levels range from 8-1 1.5 yum (black 
to blue, respectively) with a step size of 0.5 /urn. North is up and 
East is left. 

A description of observations and data reduction is given in 
Sect. [2] Models and approaches to derive synthetic observables 
are presented in Sect. [3] The stellar parameters are derived in 
Sect. [4] and compared to state-of-the-art evolutionary tracks in 
Sect. [5] A summary of the results is given in Sect. [6] 

2. Observations and Data Reduction 

Sections lZTI and 12.21 discuss the interferometric and spectro- 
scopic observations obtained with VLTI/MIDI and ISO/SWS, 
respectively. Section l2~3l discusses possible cycle-to-cycle and 
intra-cycle variability of the data. 

2.1. MIDI visibilities and spectra 

TX Psc was observed in 2004 with the 8.2 m Unit Telescopes 
and in 2005, 2006 and 201 1 with the 1.8 m Auxilia ry Telescopes 
of the Very Large Telescope Interferometer MIDI dLeinert et al.l 
l200l . MIDI covers the jV-band and provides spectrally dis- 
persed visibilities, differential phases and fluxes (resolution R = 
230 for observations in 2004, R = 30 for observations in 2005, 
2006 and 2011). 

The journal of available MIDI observations is given in Table Q] 
(electronic version only). The Mv-coverage i s plotted in F ig.Q] 
Data are reduced using M IA+EWS 1.7. lQ dJaffdl2004t Igatzka 
120051: iLeinertet all [2004). Observations are not used if one or 
more of the selection criteria discussed in iKlotz et al.1 (l2012al) 
are violated (observations that are not boldfaced in Tablefl]). 
Uniform-disk angular diameters and IRAS 12 fim flux of the cal- 
ibrator targets are given in Table |2](electronic version only). 
As only one suitable calibrator is available per observation, a 
standard multiplica tive error of 10 % is assumed for the cali- 
brated visibilities (IChesneaul 120071) . Some calibrated visibili- 
ties at baselines shorter than ~30 m are significantly noisier and 
sometimes larger than unity after 1 1.5 //m. Therefore, in the fol- 
lowing, spectro-interferometric observations are considered only 
for wavelengths shorter than 11.5 fim. 

For the calibrat ed spectra additional selection criteria are applied 
dCh esnea ul2007h : (i) the airmass difference between science and 

1 http : //www . strw . leidenuniv . nl/~ j af f e/ews/MIA+EWS- 
Manual/ index . html 



calibrator observations is < 0.2, (ii) the spectral type of the cali- 
brator is not later than M0. Considering these criteria, five spec- 
tra are derived (flagged with 'a' in Table[T). 
TXPsc is almost unresolved for baselines shorter than 32m. 
Therefore, we expect that the star is unresolved by the single- 
dish UTs. Thus, most of the mid-infrared (mid-IR) flux is located 
within the field-of-view (FoV) of the UTs and consequently also 
in the FoV of the ATs and ISO, which makes the spectra fully 
comparable to each other. 

2.2. ISO spectra 

Three spectra of TX Psc were observed with the S hort- 
Wavelength-Sp ectrometer ( SWS, Ide Graauw et al.l fl 996 ) on- 
board of ISO (ISloan et al.l 120031: (j0rgensen et al.1 l2000l> . Two 
spectra have a resolution of R ~ 200 and range from 2.36 - 
45.35 /vm. The other spectrum ranges from 2.45 - 45.20 /im and 
has a resolution of R ~ 2 000 which has been binned to the res- 
olution of the other spectra. For the ISO spectra a multiplicative 
error o f +10% is assum ed from 2.38-4.05 /im and ±5% after- 
wards (ISloan et al.ll2003l) . 

2.3. Interferometric and spectroscopic variability 

In the left panel of Fig.[2]the calibrated MIDI spectra are over- 
plotted to the ISO/SWS spectra. The flux level of the MIDI ob- 
servations is the same (within the error bars) as the ISO/SWS 
spectra taken at 1996 Nov 26 and 1997 Dec 11. This suggests 
that no significant cycle-to-cycle variation is expected in the 
mid-IR. The right plot of Fig. [2] shows part of the lightcurve 
of TXPsc in V taken from AAVSO. The spectra from 1996 
Nov 26/1997 Dec 11 were observed at a 'local' visual maxi- 
mum/minimum, respectively. In the following we will use these 
two spectra to derive the parameters of the star at different 
phases. As we do not have any information on the V magnitude 
of the spectrum from 1997 May 24, this spectrum is not used 
in the subsequent sections. Additionally, the phase for all the 
MIDI observations is unknown. However, AAVSO visual esti- 
mates from one specific observer revealed that none of the MIDI 
observations were taken at visual minima or maxima. 
Simulations with dust-free model atmospheres show that the ef- 
fect of pulsati on on the Af-band U D-radius is smaller than 0.1 
mas at 1 kpc (Paladini et al. 2009). For a star at the distance of 
TX Psc this would be ~0.3 mas. This value is beyond the limit 
of resolution of MIDI. Therefore, we do not expect any observed 
intra-cycle variability effect that is larger than the errors. This al- 
lows to combine all calibrated visibilities for the further analysis. 

3. Model description 

In the following we present the different classes of models and 
the derivation of synthetic observables that will be compared to 
observations in Sect. [4] 

3.1. Hydrostatic models 

Observed spectra and visibilities are compared to the grid 
of spherical hy drostatic model atmospheres and spectra 
of lAringer et al l (|2009). These models are computed with 
COMARCS and are generated assuming hydrostatic local ther- 
mal and chemical equilibrium. The molecular and atomic opaci- 
ties are treated in the opacity sampling (OS) approximation. The 
parameters that characterize a model are: effective temperature 
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Fig. 2. Left: MIDI flux (error bars) for different dates. Overplotted are the three ISO/SWS spectra (full lines) of TXPsc. 
AAVSO V-band light curve of TXPsc. Vertical lines mark the observation date of the ISO spectra. Colors are the same as in 
figure. 
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T e ff , metallicity Z, surface gravity g, mass M, and carbon to oxy- 
gen ratio C/O. 

For this work we limit the sample to models having solar metal- 
licity as the effect of metallicity is expected to be small for low- 
resolution spectroscopy. Additionally, there is no indication for 
largely non-solar metallicity from other properties of this star. 
The spectra cover the following parameters: 2400 < T e ff < 
4000 K with steps of 100 K; Z/Z Q = 1; -1.0 <log(g[cms" 2 ]) < 
+0.0; M/M = 1,2; C/0= 1.05, 1.10, 1.40,2.00. In order to 
get precise estimates of T e ff, additional model atmospheres 
were produced resulting in a grid spacing of AT=10K. All the 
main molecular opacities typica l for C- stars were includ ed: CO 
(iGoorvitch & Chackerianl|l994l). Co Jbuerci et al.lll974T) . HCN 



(Harris et alj|2006l) . CN ( J0rgensen|[i997l) in t he form of line 
lists, while C 2 H 2 and C 3 (IJorgensen et al.ll 19891) as OS data. CS 
is not included due to the lack of line lists and OS data. Synthetic 
spectra with a resolution of 18 000 are computed in the wave- 
length range 0.8 - 25 fim. The spectra are convolved in order to 
get the same resolution as the observed data. 
Among the output of the spherical radiative transfer code COMA 
is the monochromatic spatial intensity profile. This profile is 
used to calculate a synthetic visibility profile in the mid-IR for 
a subset of the models in the grid. A detailed des cription of the 
comp utation of the visibility profiles is given in Paladini et al. 
(120091) . 

3.2. Evolutionary tracks 

The luminosity and effective temperature that are determined 
from hydrostatic models are compared to thermally-pulsing (TP) 
AGB evolutionary tracks from Marigo et al. (in prep). We se- 
lected TP-AGB sequences with an initial scaled-solar chemical 
composition (Z = 0.014, Y = 0.273), where Z and Y denote 
the mass fractions of metals and helium, respectively. TP-AGB 
evolutionary calculations are carried out from the first thermal 
pulse - extracted fro m the PARSEC database of stellar models 
(Bres san et al.ll2012h - to the complete ejection of the envelope 
due to stellar winds. The TP-AGB tracks are based on numerical 
integrations of complete envelope models in which, for the first 
time, molecular chemistry and g as opacities are compute d on- 
the-fly with the /ESOPUS code iMarigo & Aringerll2009l) . This 
guarantees a full consistency of the envelope structure with the 



surface chemical abundances that may significantly vary due to 
the third dredge-up episodes and hot-bottom burning. The transi- 
tion from C/O < 1 to C/O > 1 is followed accurately, in particu- 
lar in the narrow range from .95 ~ C/O « 1 .05, where an abrupt 
change in the molecular che mistry and opacity is exp ected to oc- 
cur (see figures 1 1 and 16 in lMarigo & Aringerl2009l) . This point 
is particularly relevant in the context of the present work, as TX 
Psc is found to have a surface C/O slightly above unity. 



4. Stellar parameter determination 

Stellar parameters for TX Psc were determined by a number of 
authors. A summary is given in Table|3] 

In the following two different approaches will be used to de- 
termine the stellar parameters of TXPsc: (i) in Sect. 14.11 we 
use the 'classical approach' of interferometrists where spectro- 
interferometric observations are used to determine the ef- 
fective te mperature T e g ; (ii) i n Sect. 14.21 we follow the ap- 
proach of iPaladini et al.l (1201 ll Pll hereafter) where spectro- 
scopic/interferometric measurements are compared to hydro- 
static models to determine T e s and C/O ratio and to constrain 
mass and logg. 



4. 1 . Approach 1 : Geometric models 

Interferometry is a powerful tool for constraining the morphol- 
ogy and size of stars and their environments. 
A deviation from spherical symmetry for the circumstellar en- 
vironment of TX Psc was detected by s everal authors at differ- 
ent spatial scales an d wavelengths fe.g. |Cruzalebes et al. 1998; 
iRagland et all l2006t ISacuto et al.l 1201 ltat iJorissen et all l20lTl 
Hron et al in prep.). Clearly, the circumstellar environment of 
TX Psc is very complex and most likely unrelated effects shape 
the morphology in different regions, resulting in a large variety 
of structures. To study the geometry of the innermost region of 
the sta r in the Af-band, th e geometrical model fitting tool GEM- 
FIND dKlotz et alJ2012bh was used to fit the MIDI interferomet- 
ric observations of TXPsc. A spherical UD-model is able to re- 
produce the calibrated visibilities, i.e. TXPsc can be assumed to 
be spherically symmetric in the mid-IR at spatial scales probed 
by our MIDI observations. The reason why the asymmetrical 
structures detected by other works are not observed is that they 
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Table 3. Published stellar parameters of TX Psc. 
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Fig. 3. Left: Calibrated visibilities (symbols) versus baseline length for three different wavelengths. The lines represent the best- 
fitting UD-model at the given wavelength. Right: Spectrally dispersed angular diameter plus errors from the best-fitting UD-model 
(dark grey shaded area). The light grey shaded area marks the region that is omitted for the mean diameter estimation. 



were either detected at distances outside the FoV of MIDI or at 
smaller spatial scales beyond the detection limit of MIDI. This 
is supported by the differential phase measurements of MIDI, 
which do not show any deviation from zero. The calibrated visi- 
bilities are plotted together with the best-fitting UD-model in the 
left panel of Fig.[3]for three different wavelenghts. 
The right panel of Fig. [3] shows the wavelength dispersed di- 
ameter calculated with the UD-model (dark grey shaded area). 
The star appears lar g er be tween 8 and 9 //m. According to 
Gautschv-Loid leTaT] (120041) CS is affecting this wavelength 
range. Because of this molecular contamination, this region is 
omitted and a mean angular diameter of 6 = 10.51 + 0.70 mas is 
calculated by averaging the diameter from 9-1 1.5 fim. 

4.1 .1 . Effective temperature 

The temperature can be determined using the apparent bolomet- 
ric magnitude m\, \ and the angular Rosseland diameter # ross . The 
apparent bolometric magnitude m\, \ = 2.26 mag was derived by 



F. Kerschbaum by fitting a combination of blackbodies to near- 
IR and IRAS datifl 

Various de finitions for the radius can be found in literature (c.f. 
reviews by laschek et alJll991t IScholzll2003l) . where the most 
commonly used radius in atmospheric modeling is the Rosseland 
radius. It is defined by the distance between the center of the 
star and the layer having Rosseland optical depth T loss = 2. 
This radius, however, is not an observable quantity and observed 
radii have to be converted by using model considerations. In 
the following we will derive this conversion factor for hydro- 
static C-stars by using a sub set of the hydrostatic models in the 
grid of lAringeretZI(b009b to derive a mean UD-radius in the 
mid-IR (9 - 11. 5 /mi). This mean UD-radius is plotted versus 
the Rosseland radius of the corresponding hydrostatic model in 
Fig-El There is a clear correlation between the two radii, yielding 

R Ross = 0.95 Rud • (1) 



2 Method described in Kerschbaum & Hron d!996af) and reference s 
therein; near-IR data from the IRAS catalogue an d lFouaue et al. (1992). 
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rived mid-IR UD-radius of the same models (crosses). The blue 
line is a linear fit to the data. 

Table 4. Stellar parameters derived from observations (middle 
block) and calculated (right block) using approach 1 . 
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This implies that for a hydrostatic C-star the Rosseland radius 
can be approximated by the mid-IR UD-radius if the derived cor- 
rection factor of 0.95 is applied. 

Applying this correction factor yields an angular Rosseland di- 
ameter Sross = 9.99 mas. Together with the bolometric magni- 
tude the distance-independent effective temperature can be de- 
rived. Using the three different distance estimates d that are 
available for TXPsc (see Tabled a linear radius and luminosity 
can be calculated. Table|4]lists the derived and calculated stellar 
parameters of approach 1. Errors on the temperature are deter- 
mined by using the errors on #ud and by assuming an arbitrary 
error of +0.1 for mbd that accounts for the stellar variability and 
the fitting error. These values are in agreement with those given 
in literature (see Table[3). 

4.2. Approach 2: Hydrostatic models 

The short wavelength part of the ISO spectrum of TX Psc is 
dominated by the 3 /im feature which has contributions from 
HCN and C2H2. The 5 fim feature, on the other hand, is due 
to C3 and the fundamental band of CO. The region from 7 
to 8//m is dominated by H CN and C2H2 and according to 
Gautschv- Loidletal.1 (l2004h the region from 8 to 9 fim is af- 
fected by CS opacity. The low variability in the V-band (AV ~ 
0.4 mag) justifies the use of hydrostatic models in the near- as 
well as mid-IR. Additionally, hydrostatic models are able to re- 
produce large parts of the overall ISO as well as MIDI spectra. 
This indicates that the circumstellar environment of TX Psc is 



optically thin and contains almost no dust. 
In Sect. l4.2.T1 and 14.2.21 low-resolution spectroscopic observa- 
tions are compared with synthetic spectra of hydrostatic models 
to fix the fundamental stellar parameters C/O ratio and T e ff . The 
overall energy distribution as well as the bands of the molecules 
that are pre sent from 2.3-6 um put strong constrai nts on these 
parameters (ttorgensen et al-lBOOOt iLoidl et al.lfcOQlt PI 1). 
Low resolution spectroscopy does not allow to ascertain mass 
and \ogg (Fig. 6-10 in Pll). To determine these parameters 
Sect. 14.2.31 follows the approach described in Pll: spectro- 
interferometric observations are compared to models of fixed 
r e ff and C/O ratio but varying log g and mass. 

4.2.1. C/O ratio 

J0rge nsen et al.l (l2000h found the ratio between the 3 fim feature 
(HCN and C 2 H 2 ) and the 5.1 fim feature (CO and C 3 ) to be a 
sensitive measure of the C/O ratio. In order to be independent of 
distance each model spectrum is normalized to t he ISO flux at 
2.9 fim (local minimum of molecular absorption, Aringe r et al.l 
120091) . Ax 2 test is applied between 2.9 - 6.0/im to compare the 
observed ISO spectra (1996 Nov 26 and 1997 Dec 11) to the 
models. For both ISO spectra we find that the best solution is 
obtained with a C/O ratio of 1.05. Considering all solutions ly- 
ing within the 68% confidence level a C/O ratio of 1.1 can be 
defined as an upper limit. Due to the coarse grid spacing for 
C/O this value is not a strict upper limit and no lower limit can 
be assigned. Figs. [5] and [6] show the ISO/SWS spectrum plotted 
together with models of different temperatures and C/O ratios. 
The upper panels draw the region around the 3 fim and 5.1 fim 
feature to a larger scale. These plots demonstrate that a synthetic 
spectrum with a C/O ratio larger or equal 1 .4 is not able to re- 
produce the observations, because it over-evaluates the depth of 
the 5.1 fim feature. This finding is consistent with the C/O ratios 
given in literature (see Table|3]l. 

4.2.2. Effective temperature 

We expect the 3 fim feature to be a good temperature indicator 
for hydrostatic stars (PI 1). Consequently, in orderto find the best 
temperature for TX Psc, ax 2 test i s use d t0 compare the observed 
ISO spectra and the model spectra between 2.9 - 3.6/vm. Only 
model spectra lying within the confidence level of C/O (1.05, 
1.1; see Sect. 14. 2. Ti l are used for this test. 

The large grid of models allows to determine T e ff very pre- 
cisely. All solutions lying within the 68% confidence level have 
a temperature of 3 080+™ / 3170+™K for the visual mini- 
mum/maximum, respectively. The upper left panel of Figs.|5]and 
|6]show that models with higher/lower temperatures do not repro- 
duce the depth of the 3 fim feature. 

The temperature derived with approach 1 in Sect. l4.1.T1 is within 
the errors of the temperature that is derived here for the visual 
minimum and maximum. 

We confirm the finding of J0rge nsen et al.l (l2000l) that the dif- 
ference in the two ISO spectra can be explained by a tempera- 
ture change. The lightcurve suggests that this difference is due 
to variability effects, but time-series spectroscopy is needed to 
confirm this finding. 

4.2.3. Surface gravity and mass 

Interferometric observations are compared to synthetic visi- 
bilities of hydrostatic models of varying logg and mass us- 
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Fig. 5. ISO/SWS spectrum of TX Psc at visual minimum from 1997 Dec 1 1 (black line) plotted with hydrostatic models (colored 
lines) of different temperature and C/O ratio. The best fitting model is plotted in green. Model spectra are normalized to the flux of 
the corresponding ISO spectrum at 2.9 //m. 
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Wovelength [fim] 

Fig. 6. Same as Fig. [5] but for at visual maximum from 1996Nov21. 



ing a x 1 tes,; - All synthetic visibilities are computed for the 
best-fitting values of T e ff and C/O-ratio from Sects. |4~2~T1 and 
14.2.21 Synthetic visibilitie s are derived at 3 different distances: 
275 pc, 280 pc, 3 15 pc Ivan Leeuwenl l2007t IClaussen et al.l 
119871: (Bergeat & Chevallierl 120051 respectively). In Fig.|7] the 
wavelength-dispersed calibrated visibilities are plotted together 
with the synthetic visibilities of the best-fitting models. The re- 
gion between 8 - 9yum is not considered in the fitting, as the 
hydrostatic models do not include the CS opacity data. 
The middle block of Table[5] gives the best-fitting stellar pa- 
rameters for a given distance that were determined using ap- 
proach 2. It is clear from the table and from Fig.|7]fhat, given 
the error bars on the visibilities, there is a degeneracy between 
logg and mass. As the distance defines the level of visibility, 



there is also a degeneracy between distance an d logg. From 
current TP-AGB evolutionary calculations (e.g. Karakas et al. 
l2002l:lMarigo & Girardill2007l Marigo et al. in prep.) we expect 
that a 1 M TP-AGB star with solar metallicity does not make 
the transition to the C-rich domain. This suggests that the log g 
values in Table[5]found for M — 2 M are the more reliable ones. 
But, considering the degeneracy and the limited mass sampling 
(M =1,2 M ) in the grid, also models with higher masses would 
reproduce the observed visibilities. To support this statement, 
one additional model with 3 M Q is calculated and overplotted in 
Fig.|7]as dashed line. The models are almost indistinguishable. 
The right part of Table|5]gives stellar parameters that are calcu- 
lated from the derived log g, mass and effective temperature. The 
luminosities are significantly larger than the L = 5 200 L Q used 
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Fig. 7. Wavelength-dispersed calibrated visibilities plus errors (dark grey shaded area) plotted with the best-fitting synthetic visibil- 
ities of hydrostatic models (full lines). To show the degeneracy between mass and log g one model with 3 M is overplotted (dotted 
line). 



bv lGautschv-Loidl et all d2004. but comparable to L = 7 700 L e 



derived bv lClaussen et alj (11987b . The luminosities and radii are 
also in perfect agreement with the ones determined with ap- 
proach 1 (see Sect. 14. U . 



4.2.4. Photometric constraints on the best model 

The best fitting hydrostatic models of visual mini- 
mum/maximum (model with lowest Xi ogg M m Tabled are 
overplotted to the ISO and MIDI spectra as well as to pho- 
tomet ric measurements in Fig.|8] Photometric measurem ents 
from Johnson et al. il%6|), |Mendoza V. & Johnson! (1 19651) and 
ICatch nole et al.l d 19791) were observed with the Johnson filter 
system. Zero points to convert these measu rements from mag- 
nitudes to Jansky are taken from lCoxl (l2000|). These zero points 
are also used to convert obs ervations from lBergeat et al.l (11976) 
an d lBergeat & Lunell (Il980h as the authors claim that their filter 
system is similar to the Johnson filter system. Ze ro points for 
2MAS S photometry (fCutri et al.ll2003 | ) are given inlCohen et al.l 
(120031) . lOlofsson et al. dl993l) and TKerschbaum et alj d!996bf 
used th e E SQ filter syst e m. Ze r o points are tak e n fromlLe Bertrd 
d!988l) and lWamstekerl (Il98ll) . iNoguchi et af] (Il98ll) used their 
own filter system and corresponding zero points are given in 
their paper. 




Photometry: 

* Mendozo & Johnson 1965 
+ Johnson et al. 1966 
? Bengeat et al. 1976, 1980 
n Catchpole et al. 1979 



* Olofsson « 
A Kerschbai; 
o Cutri et a 



Wavelength [/i.m] 

Fig. 8. Best fitting hydrostatic models for the visual minimum 
(black line) and maximum (grey line). Superimposed are the ISO 
spectra of visual minimum (green line) and visual maximum 
(blue line), MIDI spectra (violet lines) and photometric mea- 
surements of different works (orange symbols). Model spectra 
are normalized to the flux of the corresponding ISO spectrum at 
2.9 //m. 



5. Comparison with evolutionary tracks 

We follow the approach described in PI 1 and compare the stel- 
lar parameters with new evolutionary tracks of thermally pulsing 
AGB stars (Marigo et al., in prep). 

Figure|9] depicts evolutionary tracks in the region of AGB stars 
in the H-R diagram. Overplotted are the determined luminosity 
and temperature for TXPsc for the two approaches. Only the 
best-fitting luminosity (see Table [3]l at d=280pc is plotted for 
TXPsc. Errors for the luminosity are assumed to lie in the order 
of ~40% (uncertainty on the given distance measurement). 
We note an encouraging agreement between the observed lo- 
cation of TX Psc in the H-R diagram and the predicted ranges 
of luminosity and effective temperature for a carbon-rich star 



with solar- me tallicity and C/O within a narrow interval (i.e. 
1<C/0<1.1). As we see in Fig. |9l the TP- AGB evolutionary 
tracks in the C-rich regime extend to much lower effective tem- 
peratures than the derived values for TX Psc. This cooling is 
mainly driven by the increase of the C/O ratio after each dredge- 
up episode, as well as by the progressive strengthening of the 
mass-loss efficiency. The relatively warm effective temperature 
of TX Psc suggests that this star is observed close after the tran- 
sition into the C-star domain, in an early 'quiet' stage in which 
the strong wind has not yet developed. This picture is nicely sup- 
ported by the observational findings already discussed in the pre- 
vious sections. 

It is visible from the evolutionary tracks that an AGB star with 
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Table 5. Stellar parameters derived from observations (middle block) and calculated (right block) using approach 2. 



d 

[pc] 


[K] 


T ff 
[K] 


c/o 


\ogg 


M 

[M ] 


% logg.M 


R 


[i Q ] 




275 


3080!™ 


3170!™ 


1.05 


-0.5 


1 


0.50 


295 


7019 


7876 








-u.z 


__ 






7H1 Q 


/o /O 


280 


3080!™ 


3170!™ 


1.05 


-0.5 
-0.2 


l 

2 


0.48 
0.53 


295 
295 


7019 
7019 


7876 
7876 


315 


3080!™ 


3170!™ 


1.05 


-0.6 
-0.3 


1 

2 


0.49 
0.53 


331 
331 


8836 
8836 


9915 
9915 



solar metallicity will turn into a carbon-rich AGB star only for 
masses around 2 M Q and higher. Also, the position of TX Psc 
in Fig. [9] suggests that the mass lies between 2 and 3M . This 
is not in agreement with the best fitting models having 1 M 
(Sect. 14. 23b . but in good agreement with the models having 
2.W . 
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6. Conclusion 

In this work we determined stellar parameters for TXPsc by 
comparing observations to geometric models dKlotz et al.l 
1201 2bl). state-of-the -art hydrostatic model atmospheres 
(I Aringer et alJ 120091) and evolutionary models (Marigo et 
al. in prep.). Two different approaches were used to fix the 
parameters: 

Al. Spectro-interferometric observations were used to de- 
termine a wavelength-dispersed uniform disk diameter. A 
correction factor for hydrostatic C-stars was derived from 
hydrostatic models to convert the UD diameter to the Rosseland 
radius, which was then used to determine the effective tempera- 
ture r e ff. 

A2. Spectroscopic measurements were compared to synthetic 
spectra from hydrostatic models to determine r e ff and C/O 
ratio. The mass and logg were constrained by comparing 
spectro-interferometric observations with synthetic visibility 
profiles from hydrostatic models. 

The main advantage of approach 1 is the distance-independent 
determination of T e ff. On the other hand, conversion of the 
UD-radius to the Rosseland radius and the use of the apparent 
bolometric magnitude introduces uncertainties. Approach 2 
allows to constrain not only T e ff, but also C/O, logg and M. 
One of the disadvantages of this technique is the unknown error 
that is introduced by the model. Additionally, the uncertainty 
in distance, that is needed to constrain logg and M, and the 
degeneracy between these two parameters, limits the accuracy of 
the parameter determination. This suggests that high-resolution 
spectroscopy is needed to fully discriminate between mass and 
logg. 

There is a very good agreement between the best-fitting hy- 
drostatic model atmosphere and observations (interferometry, 
spectroscopy and photometry). 

Our spectro-interferometric results are also an important tool 
to constrain and validate stellar AGB models, that are still 
subject to severe uncertainties. We found that present TP-AGB 
tracks with a detailed treatment of molecular opacities nicely 
reproduce the derived T e ff, L, C/O values for TX Psc. 
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Fig. 9. Zoom into the AGB region of the H-R diagram. Lines denote solar metallicity evolutionary tracks (Marigo et al., in prep.) 
and numbers indicate the mass on the early-AGB. Yellow/grey lines mark the region of carbon-rich AGB stars with C/O>1.0. Black 
lines mark the region of oxygen-rich AGB stars (C/0<1). For better visibility, the track of the 2 M© model is plotted with a dotted 
line. Different colored symbols refer to the luminosity and effective temperature determined in this work (for the two different 
approaches Al and A2 at visual minimum/maximum). 



Jorgensen, U. G., Almlof, J., & Siegbahn, P. E. M. 1989, ApJ, 343, 554 

J0rgensen, U. G., Hron, J., & Loidl, R. 2000, A&A, 356, 253 

Jorissen, A., Mayer, A., van Eck, S., et al. 201 1, A&A, 532, A135 

Karakas, A. I., Lattanzio, J. C, & Pols, O. R. 2002, PASA, 19, 515 

Kerschbaum, F. & Hron, J. 1996a, A&A, 308, 489 

Kerschbaum, F., Lazaro, C, & Habison, P. 1996b, A&AS, 118, 397 

Klotz, D., Sacuto, S., Kerschbaum, E, et al. 2012a, A&A, 541, A164 

Klotz, D., Sacuto, S., Paladini, C, Hron, J., & Wachter, G. 2012b, ArXiv e-prints 

Kurucz, R. 1993, Limbdarkening for 2 km/s grid (No. 13): [+0.0] to 

[-5.0]. Kurucz CD-ROM No. 17. Cambridge, Mass.: Smithsonian 

Astrophysical Observatory, 1993., 17 
Lambert, D. L., Gustafsson, B., Eriksson, K, & Hinkle, K. H. 1986, ApJS, 62, 

373 

Lasker, B. M., Bracker, S. B., & Kunkel, W. E. 1973, PASP, 85, 109 
Le Bertre, T. 1988, A&A, 190, 79 

Leinert, C, Graser, U., Przygodda, E, et al. 2003, Ap&SS, 286, 73 
Leinert, C, van Boekel. R., Waters, L. B. F. M., et al. 2004, A&A, 423, 537 
Loidl, R., Lancon, A., & J0rgensen, U. G. 2001, A&A, 371, 1065 
Marigo, P. & Aringer, B. 2009, A&A, 508, 1539 
Marigo, P. & Girardi, L. 2007, A&A, 469, 239 

Marti'- Vidal, I., Marcaide, J. M., Quirrenbach, A., et al. 2011, A&A, 529, A115 

Mendoza V., E. E. & Johnson, H. L. 1965, ApJ, 141, 161 

Neilson, H. R. & Lester, J. B. 2008, A&A, 490, 807 

Noguchi, K, Kawara, K, Kobayashi, Y., et al. 1981, PASJ, 33, 373 

Ohnaka, K, Tsuji, T., & Aoki, W. 2000, A&A, 353, 528 

Olofsson, H, Eriksson, K., Gustafsson, B., & Carlstroem, U. 1993, ApJS, 87, 
305 

Paladini, C, Aringer, B., Hron, J., et al. 2009, A&A, 501, 1073 

Paladini, C, van Belle, G. T., Aringer, B., et al. 201 1, A&A, 533, A27 (PI 1) 

Querci, F, Querci, M., & Tsuji, T. 1974, A&A, 31, 265 

Quirrenbach, A., Mozurkewich, D., Hummel, C. A., Buscher, D. F, & 

Armstrong, J. T. 1994, A&A, 285, 541 
Ragland, S., Traub, W. A., Berger, J.-R, et al. 2006, ApJ, 652, 650 
Ratzka, T. 2005, PhD thesis, Max-Planck-Institute for Astronomy, Konigstuhl 

17, 69117 Heidelberg, Germany 



Richichi, A., Chandrasekhar, T., Lisi, F, et al. 1995. A&A, 301, 439 

Sacuto, S., Aringer, B., Hron, J., et al. 201 la, A&A, 525, A42 

Sacuto, S., Jorissen, A., Cruzalebes, P., et al. 2011b, in Astronomical Society 

of the Pacific Conference Series, Vol. 445, Why Galaxies Care about AGB 

Stars II: Shining Examples and Common Inhabitants, ed. F. Kerschbaum, 

T. Lebzelter, & R. F. Wing, 171 
Samus, N. N., Durlevich, O. V., & et al. 2009, VizieR Online Data Catalog, 1, 

2025 

Scholz, M. 2003, in Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 4838, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, ed. W. A. Traub, 163-171 

Sloan, G. C, Kraemer, K. E., Price, S. D., & Shipman, R. F. 2003, ApJS, 147, 
379 

van Leeuwen, F. 2007, A&A, 474, 653 
Wamsteker, W. 1981, A&A, 97, 329 

Wittkowski, M., Boboltz, D. A., Driebe, T., et al. 2008, A&A, 479, L21 
Wittkowski, M., Boboltz, D. A., Ireland, M., et al. 2011, A&A, 532, L7 
Wittkowski, M., Hummel, C. A., Johnston, K. J., et al. 2001, A&A, 377, 981 



9 



Klotz et al.: Constraining stellar parameters for TXPsc, Online Material p 1 



Table 2. Properties of the calibrator targets. 



HD 


Name 


Sp. T. a 


Fn a 


7T. 








[Jy] 


[mas] 


HD 48915 


Sirius 


Al 


143.1+3 


6.08+0.03 


HD 20720 


t04 Eri 


M3/M4 


162.7±6 


10.14+0.04 


HD 224935 


YYPsc 


M3 


86.9±5 


7.25+0.03 


HD 49161 


17Mon 


K4 


10.4±5 


2.44+0.01 


HD 18884 


a Cet 


M1.5 


234.7±3 


12.28+0.05 


HD 45348 


Canopus 


F0 


154.8+3 


6.87+0.03 



Notes. (a) http : //www . eso . org/ observing/ df o/quali ty/MIDI/ qc/ 

calibrators_obs . html 

(b) http : // simbad . u-strasbg . f r/ simbad/ 
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Table 1. Journal of MIDI observations of TX Psc sorted for projected baseline length. 



# 


Object 


Date/Time 


Configuration 


Proj.Base 
[m] 


Proj. Angle 
[°] 


Obs. seeing 
["] 


Airmass 


1 


TXPSC 


2006-08-28 02:50 


E0-G0 


9.4 


81 


1.41 


1.718 


2 


TXPSC 

HD48915 


2006-10-19 00:19 
2006-10-19 08:26 


E0-G0 


11.9 


70 


0.57 


1.387 
1.325 


3 


TXPSC 

HD48915 


2006-10-18 05:46 
2006-10-18 06:57 


E0-G0 


13.2 


66 


0.89 
0.87 


1.659 
1.278 


4 


TXPSC 


2005-06-29 08:55 


E0-G0 


13.7 


82 


0.54 


1.238 


5 


TXPSC 


2010-09-07 04:13 


E0-G0 


14.4 


73 


0.93 


1.194 


6 


TXPSC 


2006-08-28 05:48 


E0-G0 


15.6 


73 


1.30 


1.136 


7 


TXPSC 


2009-11-15 02:18 


E0-G0 


15.7 


72 


1.26 


1.225 


8 


TXPSC 


2006-10-18 03:18 


E0-G0 


16.0 


73 


1.00 


1.149 




HD48915 


2006-10-18 06:57 








0.87 


1.278 


9 


TXPSC 


2006-09-20 02:28 


D0-G0 


25.1 


71 


1.21 


1.327 


10 


TXPSC 


2006-09-20 03:23 


D0-G0 


29.0 


73 


1.90 


1.189 


11 


TXPSC 


2006-10-16 02:12 


H0-G0 


30.5 


73 


0.67 


1.150 


12 


TXPSC 


2009-11-16 02:13 


HO-GO 


31.0 


72 


0.95 


1.219 


13 


TXPSC 

HD20720 


2006-09-20 05:47 
2006-09-20 09:48 


DO-GO 


31.6 


72 


1.64 
0.89 


1.201 
1.095 


14 


TXPSC 


2009-11-16 01:20 


HO-GO 


32.0 


73 


1.36 


1.146 


15 


TXPSC 


2006-09-21 01:48 


KO-GO 


44.0 


68 


1.01 


1.489 


16 


TXPSC 

HD48915 


2006-09-21 03:32 
2006-09-21 08:36 


KO-GO 


59.3 


73 


1.02 
1.51 


1.169 
1.309 


17 


TXPSC 


2006-09-17 04:20 


AO-GO 


62.0 


73 


1.36 


1.140 


18 


TXPSC 

HD48915 


2006-09-21 05:55 
2006-09-21 08:36 


KO-GO 


62.6 


72 


1.39 
1.51 


1.225 
1.309 


19 


TXPSC 


2008-09-27 04:45 


G1-D0 


63.6 


133 


1.00 


1.155 


20 


TXPSC 

HD224935 


2006-08-17 07:00 
2006-08-17 07:23 


AO-GO 


63.8 


73 


0.39 
0.43 


1.137 
1.060 


21 


TXPSC 

HD49161 


2004-10-30 03:28 
2004-10-30 09:19 


UT2-UT4 


84.3 


80 


0.78 
0.72 


1.244 
1.196 


22 


TXPSC 

HD18884 


2006-08-16 07:06 
2006-08-16 07:57 


A0-G1 


87.1 


113 


0.85 
0.81 


1.137 
1.353 


23 


TXPSC 


2006-08-15 06:20 


A0-G1 


90.3 


113 


0.88 


1.145 


24 


TXPSC 


2011-09-21 05:39 


Il-Al 


104.0 


84 


1.03 


1.186 


25 


TXPSC 


2011-10-05 04:21 


K0-A1 


129.0 


67 


0.57 


1.158 


26 


TXPSC 

HD45348 


2011-10-02 05:44 
2011-10-02 06:42 


A1-J3 


140.0 


46 


1.02 
0.88 


1.305 
1.512 



Notes. Calibrators that are used to calibrate the data are given below the science target. The baseline configuration, projected baseline length and 
position angle of the observations are given. The observatory seeing and airmass give the average atmospheric conditions present at the observatory 
during the observation. Observations that are used in the astrophysical interpretation are marked in boldface. 
(a) Calibrated spectrum available 



